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Target nucleotide sequence T 



GTAGTCAGGCCAT 



Complementary GTCAGGCCAT T 
sequence Q 



Common subsequence? 



All Hybridization 



Probe P ATGGCCTGAC 



(a) 



Subsequence S ATGCGCG-TAAGT 



I 1 



Complementary ATG- - CGTTAACT 
sequence Q 

Twice insertion 

Gnce insertion 

Once substitution 

Resemblance r = 0.6364 
Maximum edit distance k = 4 



(b) 



Figure 2 
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Figure 3 
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Receive target nucleotide sequence Ti and 
nucleotide sequence of probe P 



Determine complementary sequence and 
receive edit distance k 



-S10 



S12 



Evaluate binding possibility, 

(largest.prefix^length.with. 
common.subsequence) 



Store evaluation result and 
determine binding possibility 



,S14 



S16 



Figure 4 
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Receive target nucleotide sequences 
Tl to Tn and nucleotide sequence of 
probe P 



Determine complementary sequence, and 
receive edit distance k, chain length of 
partial sequence, and chain length of overlap 
portion . 



Set counter value t to 1 



Store evaluation of binding possibility 
(largest_prefix - length - with_common^ 
subsequence; check.target..) 



Determine evaluation result 



t=t+i 
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Read out target nucleotide sequence 



S36 



no 






.S34 












Determine binding possibility from 
evaluation result and provide user with it 



Figure 
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Figure 6 
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Initialization step 




yes 



Figure 7 
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Figure 8 



8/18 




Figure 9 
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//function largest_prefix_with_common_subsequence 

largest^prefix.with^common^subsequencelsequence T[1..n], sequence Q[1..m], k) { 
for (i=0;i<=k;i++) positions[i] =0; 
min_k = 0; 

for (probe_position=1; probe_position<=m; probe j)osition++) { 
max_k = min {k, probe_position); 
for (i=max_k; i>=min_k ;i — ) { 
if (i==0) { 

positions[i] = next_char_position(T, positions[i]+1 , 
Q[probe_position]); 
] else { 

positions[i] = min { next_char_position(T, positions[i]+1 . 
Q[probej>osition]), 

positions[i-1] ) 

) 

if (positions[i] > n) { 
min k = i+1; 

) 

} 

if (min.k > k) ( 

(J) Absence of common subsequence with length of m-k or longer 
return m+1; 

] 

1 

© Presence of common subsequence with length of m-min_k or longer 
return positions(k) 
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//function check_target 

check_target(sequence T[1..n], sequence Q[1..m],k, overlapjength, checkjength) { 
positional; 
while(position n) { 

length = largesLprefix^length^with.common^subsequenceC 

T[position..n], Q[1..m], k); 
if (length < checkjength) ( 
if (check_exactly( 

T[position..(min(n. position+checkjength))], Q. k)) ( 
position = position + checkjength; 
} else { 

return ("T binds to P") 

) 

} else if (length -1 <= n - position + 1)( 

position = position + length - overlapjength - 1 ; 

} else ( 

position = n+1; 

) 

i 

return ("T does not bind to P*'); 
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Figure 12 
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Edit distance k— 3: positionsfO] 

m 123 4. 5678 



Complementary sequence 



Target nucleotide sequence 

i 




CGCGCATGAA 
1 2 3 4 5 678 9 10 



[positions[0],m]= (positions[0], 2, 3, 5, - - - - -) 



(a) 



Edit distance k-3: positions[1] 



m 



1 2 3 4 5 6 7 8 



Complementary sequence 



Target nucleotide sequence 




1 23456789 10 



[positions [1],m]= (positions[1], 1, 3, 5, 6, 7, 8) 



(b) 
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Figure 17 
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Comparative example 
(1108.33s) 




0 1 — 1 1 1 — i 1 

70 75 80 85 90 95 100 

Resemblance (%) 
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